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Abstract. The restricted numerical range Wr{A) of an operator A acting on a D- 
dimensional Hilbcrt space is defined as a set of all possible expectation values of this 
operator among pure states which belong to a certain subset R of the of set of pure 
quantum states of dimension D. One considers for instance the set of real states, or 
in the case of composite spaces, the set of product states and the set of maximally 
entangled states. Combining the operator theory with a probabilistic approach we 
introduce the restricted numerical shadow of A - a normalized probability distribution 
on the complex plane supported in Wr{A). Its value at point z e (D is equal to the 
probability that the expectation value (^|j4|^) is equal to z, where represents a 
random quantum state in subset R distributed according to the natural measure on 
this set, induced by the unitarily invariant Fubini-Study measure. Studying restricted 
shadows of operators of a fixed size D = NaNb we analyse the geometry of sets of 
separable and maximally entangled states of the Na x Nb composite quantum system. 
Investigating trajectories formed by evolving quantum states projected into the plane 
of the shadow we study the dynamics of quantum entanglement. A similar analysis 
extended for operators on D = 2 3 dimensional Hilbert space allows us to investigate 
the structure of the orbits of GHZ and W quantum states of a three-qubit system. 
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1. Introduction 

Recent studies on quantum entanglement, a crucial resource in the theory of quantum 
information processing, contributed a lot to our understanding of this deeply non- 
classical phenomenon (see e.g. [T] and references therein). In particular, some progress 
has been achieved in elucidating the geometry of quantum entanglement [21 El HI G3 [3 E] , 
but several questions concerning this topic remain still open [HI [9j [10]. It was 
also suggested that a geometric approach is useful to classify and quantify quantum 
entanglement [H [HI [12] . 

The phenomenon of quantum entanglement, non-classical correlations between 
individual subsystems, may arise in composite physical systems. Consider then the 
simplest composite system, which consists of two parts and can be described in a Hilbert 
space with a tensor product structure, = Ha ®Hb = 7~t,N A ®Hn b - Any product 
pure state, \ip) = \4>a) ® \4>b), is called separable, while all other pure states are called 
entangled. The set S of all separable pure states has the structure of the Cartesian 
product of the complex projective spaces [8], S = CP^- 4 " 1 x CP^ 3-1 . 

Among entangled pure states of a bipartite system one distinguishes the set £ of 
maximally entangled states, such that the partial trace of the corresponding projector 
is proportional to the identity matrix. The set £ contains the generalized Bell state, 
\ip+) = Ylii Im)> an d all states obtained from it by a local unitary transformation, 
Ua <S) Ub\4>+) £ £■ The set of maximally entangled states is thus equivalent to 
U(iV)/U(l). The structure of the set £ and other sets of locally equivalent entangled 
states of a bipartite system was studied in [13]. 

In this work we propose to analyse the geometry of the set of entangled and 
separable quantum states using the algebraic concepts of the numerical range and 
numerical shadow of an operator. For any operator A acting on the complex Hilbert 
space Hp one defines its numerical range [HI [15] as a subset of the complex plane which 
contains expectation values of A among arbitrary normalized pure states, 

W(A) = {z:z = ty\A\1>), IV>> e U D , (V#> = 1}. (1) 

Due to the classical Toeplitz-Hausdorff theorem the set W(A) is convex - see e.g. [To] . 
The differential topology and projection aspects of numerical range were investigated in 

PUtH] 

The standard notion of numerical range, often used in the theory of quantum 
information [T51 I2"0j |2"T] . can be generalized in several ways [221 1231 121] • For instance, 
for an operator acting on a composite Hilbert space one defines the product numerical 
range [2S] (also called local numerical range [21]) and a more general class of numerical 
ranges restricted to a specific class of states [2T] , 

W R (A) = {z:z = (1>\A\1>), meRc tt D , = l}. (2) 

Here R C Qd denotes a selected subset of the set of pure quantum states of a given size 
D. For instance, one may consider the set of all real states, or, in the case of composite 
spaces, the set of complex product states or the set of real maximally entangled states. 
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For any operator A acting on Hd one defines a probability density Pa{z) on the 
complex plane (261 l2Zj, supported in the numerical range W(A), 



Here denotes the unique unitarily invariant measure on the set Qjj of D- 

dimensional quantum pure states, also called Fubini-Study (FS) measure. In other 
words the shadow of A at a given point z e C characterizes the likelihood that the 
expectation of A among a random pure states is equal to z. Sometimes it is convenient 
to treat the numerical shadow as a probability measure ^, on a complex plane. For 
any measurable set E it reads Ha{E) = : (V'l-^IV') £ 

If the operator A is normal, AA* = A^A, its shadow Pa{z) can be interpreted as a 
projection of the set of classical states - the D-dimensional regular simplex of probability 
distributions - into a two-plane (28]. In the more general case of a non-normal A its 
shadow can be associated with a probability distribution obtained by projecting the 
set Qd of quantum pure states of size D into a plane. Thus choosing a matrix A to 
be analysed one fixes the relative position of the set Qd and determines the direction, 
along which it is projection on the plane. Hence investigating all possible shadows of 
various operators of a given size D one gathers information about the structure of the 
set fin- 

Combining the ideas of the restricted numerical range ^ with the numerical shadow 
^ one is led in a natural way to the definition of a restricted numerical shadow, 



where R is the selected subset of £lr>. In particular, choosing the appropriate subsets 
of the set of pure states we define the numerical shadow restricted to real states or 
SU(2)-coherent states. In the case of a composite Hilbert space one defines the shadow 
with respect to separable states or maximally entangled states. The restrictions can be 
combined so one can consider the shadow restricted e.g. to real separable states, or 
in the case of operators acting on 8 dimensional space, which describes a three-qubit 
system, one may study the shadow with respect to real GHZ states or complex In- 
states. It will be convenient to use simplified terms, so for brevity we will slightly abuse 
the notation and write about 'separable shadow', 'real entangled shadow' or 'complex 
GHZ shadow' of a given operator. Using the notion of probability measure we will 
denote this restricted by \x\. 

On one hand, for a given matrix A one may study its restricted numerical shadows 
Pa{ z ) determined by a given set R of quantum states. Alternatively, investigating 
numerical shadows of various matrices of a fixed dimension D with respect to a concrete 
subset R G £Id one may analyse its geometry. In this work we study in this way the 
structure of the set of real, separable and maximally entangled quantum states for a 
two and three-qubit system. 
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This paper is organized as follows. In section 2 some basic properties of the 
numerical range and the numerical shadow are reviewed. In section 3 we discuss 
numerical shadow restricted to real states. In the case of operators acting on the Hilbert 
space with a tensor product structure, H = Hn <8> "Hat, one defines classes of separable 
and maximally entangled states. Numerical ranges restricted to these sets are analysed 
in section 4 and 5 respectively. In section 6 we show that investigations of trajectories 
formed by evolving quantum states projected into the plane of the shadow contributes 
to the understanding the dynamics of quantum entanglement. In section 7 we discuss 
the simplest case of composite system consisting of three parts and described in the 
Hilbert space "H 8 = 7-L 2 ® "%2 ® %2- m the set fl 8 of pure states acting on this space one 
defines two classes of entangled states called GHZ and W. For any operator acting on 
H$ it is then natural to introduce the shadow restricted to states locally equivalent to 
GHZ or W , and these are investigated in section 6. 

2. Standard numerical shadow and the geometry of quantum states 

2.1. Classical and quantum states 

Elements of a Hilbert space are used as basic objects of the quantum theory. A 
physical system with D distinguishable states can be described by an element \ip) of the 
complex Hilbert space %d- It is assumed that such a pure quantum state is normalized, 
\\4>\\ 2 = (i^l^) = I? so it belongs to the hypersphere of dimension 2D — 1. One identifies 
any two states, which differ only by a global phase, ~ \<fi) = e ia \ip). The set of all 
pure quantum states Qd acting on %d is therefore equivalent to the complex projective 
space, Q D = CP D_1 (see e.g. [8]). 

Any convex combination of projectors onto pure states forms a mixed quantum 
state, p = X/iPilV'iKV'il with Pi > and YliiVi — 1- Any such state p, also called a 
density operator, is Hermitian, positive and normalized by the trace condition, Trp = 1. 
Thus the set Qd of all density matrices of order D forms a convex set of D 2 — 1 real 
dimensions. The projectors corresponding to the extremal points of the set of density 
operators form the set of pure quantum states, Qd = CP D_1 C Qd- In the one qubit 
case, D = 2, the set of pure quantum states forms the Block sphere, VL2 = CP 1 , which 
in this case is equivalent to the boundary of the Bloch ball Qi. For a larger dimension 
the {2D — 2) -dimensional set of pure states Qd forms only a zero-measure subset of the 
D 2 — 2 dimensional boundary of the set Qd- 

Note that the definitions of the set of pure and mixed quantum states are unitarily 
invariant, so they can be formulated without specifying any basis in the Hilbert space. 
On the other hand, among quantum states, one distinguishes also the set of classical 
states, which are formed by the density matrices diagonal in a certain basis. Hence any 
classical state p is represented by a normalized probability vector, p = {xi, x 2 , ■ ■ ■ , xd) 
such that Xi > and Ylf=i Xi — 1. The set of classical states forms thus the probability 
simplex Ad-i C R d_1 . Each of D corners of the simplex represents a classical pure 
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state, and their convex hull is equivalent to Ad_j. In the one qubit case, D = 2, the set 
Ai of classical states forms an interval which joins two poles of the Bloch sphere and 
traverses across the interior of the Bloch ball of one-qubit mixed quantum states. 

2.2. Numerical range and the set of quantum states 

The definition Q of the numerical range of a matrix A of order D specifies W (A) as the 
set of expectation values (ip \A\ip) for normalized pure states In fact the compact set 
W(v4) in the complex plane can be considered as a projection of the set Qc of quantum 
states onto a two plane. The following facts were recently established in [28] . 

Proposition 1 Consider the set of classical states of size D, which forms the regular 
simplex A^_i in R D_1 . Then for each normal matrix A of dimension D there exists an 
affine rank 2 projection P of the set Ajj^i whose image is congruent to the numerical 
range W(A) of the matrix A. Conversely for each rank 2 projection P there exists a 
normal matrix A whose numerical range W(A) is congruent to image of Ao-i under 
projection P. 

Proposition 2 Let Qd denote the set of quantum states of dimension D embedded in 
H D ~ l with respect to Euclidean geometry induced by Hilbert- Schmidt distance. Then 
for each (arbitrary) matrix A of dimension D there exists an affine rank 2 projection P 
of the set Qo whose image is congruent to the numerical range W(A) of the matrix A. 
Conversely for each rank 2 projection P there exists a matrix A whose numerical range 
W(A) is congruent to image of Qd under projection P. 

Thus by fixing the dimension D and selecting various operators A of this order and 
analysing their numerical ranges W (A) we may gather information about the projections 
of the set of quantum states. For instance, for D = 2 the projections of the Bloch sphere 
form ellipses (which could form circles or may reduce an intervals), and so the numerical 
range of any operator A of size D = 2 may be the result. Similarly the numerical range 
of operators of size D = 3 can be viewed as projections of the set Sl 3 = CP 2 of single 
qutrit quantum states, while the numerical range for matrices of order D = 4 can be 
associated with projections of the set ^4 = CP 3 [26l 128] . 

2. 3. Measure on the space of quantum states and the numerical shadow 

The Haar measure on the unitary group U{D) induces on the set Qd of quantum pure 
states the unitarily invariant Fubini-Study measure fi. In the simplest case of D = 2 this 
measure corresponds to the uniform distribution of points on the Bloch sphere Q2 = S 2 . 

For any operator A of size D we may compute its expectation value for random pure 
state e Ofl chosen with respect to measure fi. In this way for any operator A we 
define its numerical shadow - a probability distribution ^ on the complex plane. Note 
that by construction the numerical shadow Pa{z) is supported on the numerical range 
Furthermore, the numerical shadow is unitarily invariant: Pa{z) = Pjjau^{ z )- 
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For any normal matrix A, which commutes with its adjoint, the numerical shadow 
Pa{ z ) covers the numerical range W(A) with the probability corresponding to a 
projection of a regular simplex Ad-i of classical states embedded in R D_1 ) onto a 
plane. In general, for a non-normal operator A acting on Hd, its shadow covers the 
numerical range with the probability corresponding to an orthogonal projection 

of the complex projective manifold Qd = CP D_1 onto a plane. 

In this work we will analyse the numerical range restricted [21] to a certain subsets 
R of pure quantum states and the corresponding restricted numerical shadows, with 
probability density determined by random states distributed on the subset R according 
to the measure induced by the Fubini-Study measure on fl D . 

2.4- Standard numerical shadow for a diagonal matrix 

Consider a diagonal matrix of size D, namely X = diag(xi, . . . xd). Let |1) be an 
arbitrary fixed pure state in Hd. Then the entire set of random pure states can be 
obtained as = U\l), where U G U(D) is a random unitary matrix distributed 
according the Haar measure. Thus the expansion coefficients of the state |?/>) read 
(U U ,U 12 ,...,U 1D ). 

The numerical shadow of the operator X is defined as the density distribution of 
random numbers z := (ip\X\ifj), where is a random state defined by the random 
unitary matrix U. Therefore 

D 

z:= (l\U*XU\l) =Y^x l \U ll \ 2 = x-q (5) 

Note that the vector of coefficients q = qi, . . .qo with q { = \Uij\ 2 belongs to the 
D — 1 dimensional simplex A^-i- So the complex random variable z can be treated as 
a scalar product, z = x ■ q where the complex vector x is given by the diagonal of the 
operator X, while the probability vector q is random. Its distribution inside the simplex 
depends on the distribution in the space of unitary matrices U(-D). In cases of interest 
for us this distribution belongs to the class of Dirichlet distributions parametrized by a 
real number s, 

D 

P s (q) oc *(5>-l) (Qiq2---qD) s - 1 ■ (6) 

i=i 

In particular, we are interested in two situations. 

A) Standard (complex) numerical shadow j26j [28], generated by the Haar measure on 
U(-D). Then the vector q is distributed uniformly with respect to the flat (Lebesgue) 
measure on Ad-i, i.e. the Dirichlet measure with the Dirichlet parameter s — 1 

®- 

B) Real numerical shadow, generated by the Haar measure on the group O(N) of 
orthogonal matrices. Then the vector p is distributed with respect to the statistical 
measure, Pi/2{q) oc l/s/qi ■ ■ ■ qD, i-e. the Dirichlet measure with the Dirichlet 
parameter s = 1/2. 
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In both cases the shadow covers the entire numerical range of X, equal to the convex 
hull of the spectrum {x{\f =l with a density determined by the appropriate Dirichlet 
distribution. The case of real numerical shadow, i.e. the shadow with respect to the 
real pure states is treated in more detail in the subsequent section. 



3. Numerical shadow with respect to real states 

Quantum states belonging to a complex Hilbert space form a standard tool of the 
quantum theory. However, in some cases it is instructive to restrict attention to real 
states only. On one hand the set of the real states is easier to analyse than the full set 
of complex states of larger dimensionality. For instance, the phenomenon of quantum 
entanglement can be studied for the case of real states of a four-level systems, sometimes 
referred to as a pair of rebits (i.e. real bits) [29] . 

On the other hand is some physical applications it is easier to use real orthogonal 
rotation matrices to construct elements of the entire set of real density matrices. 
Therefore in this section we will study the shadow of an operator with respect to the 
set of quantum states with real coefficients. Such an investigation allows us to improve 
understanding of the structure of the set of real mixed states and its natural subsets. 



3.1. Real shadow of operators of size D 



2,3,4 



Consider the subset of pure quantum states of size D, that can be represented by 
real expansion coefficients in a given basis. This set forms a (D — 1) dimensional real 
projective space RP^ 1 C CP D_1 . 

Substituting into the definition ^ of a restricted numerical range [TS], EI] the set 
RP D_1 C CP D_1 as the subset R of Qd we arrive at the numerical range restricted to 
real states. In general the restricted numerical range needs not to be convex or simply 
connected. Similarly, restricting the integration in Q to the set RP D_1 of real quantum 
states we obtain the shadow of the operator A with respect to real states. For brevity 
we will also use a shorter expression, the real shadow of A. 

In the case of D = 2 the set of real states forms the real projective space equivalent 
to a circle, RP 1 ~ S 1 . Thus the real shadow of a generic, non-normal matrix of size two 
forms a singular probability distribution supported on an elliptic curve in the complex 
plane, and not inside its interior. 

The real shadows presented in Fig. [T] are obtained for illustrative operators of size 
D = 2 
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(9) 



These probability distributions can be thus interpreted as shadows of real projective 
spaces HP 2 , HP 3 and HP 4 on a plane. In the case 1(b) and 1(d) the shadow is supported 
on a real line, so we plot the corresponding probability distribution Pjf (x). In the other 
cases 1(a), 1(c) and |l(e) the real shadow is supported on the complex plane, so the 



density is encoded in the grey scale. The circle drawn by red dotted line represents 
the image of the sphere of dimension D 2 — 1, in which the set of quantum states can 
be inscribed. The blue dotted line represents the boundary of the standard numerical 
range of an operator. The restricted shadow is supported in this set or its subset. 

The numerical shadow carries also some information about the higher rank 
numerical range [301 EI] of an operator. For instance, the darker area of the shadow 
corresponding to a larger probability [28] allows one to recognize in Fig. [2] the numerical 
range of rank 2, written A 2 (U) of selected unitary matrices. In the case of D = 4 shown 
in Fig. [2| |(a)| it is equal to a single point at which the two diagonals of the quadrangle 
cross, while for the case D = 5 shown in panel (b) the numerical range of rank two 
is represented by the inner pentagon located inside the numerical range in this case 
numerical shadow appears to be uniform in the set A2(t/s). In the case D = 7 shown 
in the panel 



one can see darker areas between segments connecting every other 



eigenvalue. 



3.2. A general approach to real shadow 

Given M e Mjv(C), the real shadow is taken to mean the distribution of (Mx,x) 
when x e R N is a unit vector randomly chosen from the uniform distribution on 
S"- 1 = {xeR N : \\x\\ = 1}. 

If M itself is real, we have (Ms, x) = (x, Mx) so that the distribution of (Mi, x) is 
the same as that of the symmetric (M + M T )/2. Since the uniform distribution on S N_1 
is invariant under orthogonal transformations, we are free to diagonalize this symmetric 
matrix. In other words, the real shadow of M e Mn(M) is that of diag(ai,. . . , a/v), 
where are the (real) eigenvalues of (M + M T )/2. Thus we need the distribution of 

N 

E a ^- ( 1Q ) 
k=i 

It is known that x uniform on S N ~ l implies that p = (xf, . . . , a; 2 ) has a Dirichlet- 
(1/2) distribution on the simplex Aat_i = {p E [0, 1]^ : J2i Pk — 1}> its density is 
proportional to 1/ ■Jp^pi ■ ■ ■ Pn- 
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Figure 1. Numerical shadow with respect to real states of illustrative operator of size 
D = 2 — panel (a) represents projection onto the plane of real projective space HP 2 , 
D = 3 — panels (b) and (c) represent projections onto the plane of real projective 
space HP 3 , while the real shadows of operators of D = 4, panels (d) and (e) represent 
projections of real projective space HP 4 . Dashed circles represent the image of the 
outspherc in which the set of quantum states is inscribed while dotted lines denote 



the standard numerical range. Plots (a) (c) and (e) are done for matrices translated 



in such a way that their traces are equal to zero and suitably rescaled as described in 

EHl. 



The distribution of (10) is known in general [36J, so for < Ojv-i < • • ■ < fli we 

have 

As an example, let us consider the Jordan nilpotent J3 = ( 00 1 ^ . For the real 

shadow of J3, we first note that it is the same as that of diag(l/v / 2, 0,-l/V2). Rescaling 
for convenience, we compute the real shadow of diag(l, 0, — 1). The density gs{t) of 
Pi — p 3 is outside [—1, 1], satisfies g?,(—t) = gs(t), and, in view of the Dirichlet-(l/2) 
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(a) diag({ e 2 - ifc / 4 } 4 fc=1 ) (b) diag({e 2 - fc / 5 }^ =1 ) (c) diag ({e 8 "*/ 7 }^) 



Figure 2. Real shadows of diagonal unitary matrices U of size D having eigenvalues 
evenly distributed on the unit circle, |(a)| D = 4, |(b)| D = 5 and (c) D = 7. The 
shadow is supported in the regular polygon - the numerical range Ax(U) = W(U). 
Inner polygons represent higher order numerical ranges A2(U) and As(U). 



distribution on A 2 , g^it) is proportional to 

* , t (12) 

y/(x + t)(l - t - 2x)x 

for < t < 1. In particular, #3(0) = +00. In |36j the case iV = 3 is treated in terms of 
hypergeometric functions and we can compute 

93(i) = 27S 2Fl (^ ;1 '^r)' (13) 

for t E (0,1). 

For a general M e M N (C), one may write M = A + iB with A, B e M N (R), and by 
diagonalizing (A+A T )/2 and (B+B T )/2 we may understand the distributions of (Ax, x) 
and (Bx, x) separately, using the techniques above. Usually the joint distribution will 
be obscure, but we may have some hope of understanding it in certain cases. If, for 
example, M is normal and presented in diagonal form, it may be possible to see how 
the horizontal and vertical distributions knit together. 

Fig. [3] shows various shadows of M when M is a 4 x 4 unitary (in diagonal form). 



The sort of ID shadow density seen in the left panel of Fig. 3(a) seems typical; notably, 
the real density of a 4 x 4 Hermitian is constant between the two middle eigenvalues. 
Presumably, the vertical density has a similar form. Fig. [4] presents two examples of 
real shadow of generic unitary matrices of dimension five. Support of the shadow does 
not coincide with numerical range of those matrices because a generic unitary matrix 
can not be diagonalized using only orthogonal matrices. 
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(a) Complex shadow (left plot) and real shadow (right plot) of 
a Hermitian matrix diag(0, 1, 3, 5). 




Figure 3. Comparison of real and complex numerical shadow for Hermitian (panel 
3(a) I and unitary (panel [3(b)] ) matrices. 
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Figure 4. Real shadow of two generic unitary matrices of dimension five is supported 
on the subset of the numerical range formed by the pentagon plotted. 



3.3. Real shadows and their moments 

Since methods based on moments of the shadow distribution were so effective in the 
complex case, we may try to mimic them for the real shadow. It is not hard to adapt 
the method outlined at the beginning of the proof of Proposition 5.1 from to show 
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that, for a E 



R(a) 



x a = 



sh- 



unless each component of a is even, and that 



N /i* 



i?(2/3) 



n£=i 



•T)[j9| 



(14) 



Here we use the shifted factorial notation: (x) n = x(x + 1) . . . [x + n — 1), with the 
convention (x) = 1; also \(3\ denotes £)i Ak- 
in principle, we can use Eq. ( 14 ) to evaluate the moments of the real shadow density 
g(t) of diag(ai, ...,a N ) 



J t n g(t) 



dt 



S N-1 



akX k = Yl 



(15) 



In the complex case we were able to progress beyond this point, obtaining such effective 
relations for the moments as the determinant relation (see: Eq. (10) in [26\). 

As a test case, we may try to use Eq. (15) to find the moments of gs(t), the real 
shadow density of diag(l, 0,-1), discussed in the previous example. We have 

f t n 93 (t)dt= f (xi-xir=j2(t) { - 1)kR{2{n - k)A2k) 

J-x Js 2 k=0 \ J 



and finally we obtain that 



/! n 
t n g,(t)dt = Y, 
1 n 



(16) 



It can be shown, that for even n = 2m we have 



(2m + l)(2m + 3) . . . (4m + 1) ' 



To see this we start with equation Eq. (16). Changing the index of summation from k 
to n — k shows that the sum equals (— l) n times itself, hence equals zero for odd n. Now 
suppose that n = 2m, m = 0, 1, 2, . . .. The sum can be written in hypergeometric form 



1 (-) 
t n gz(t)dt X ' 2 ' 



l\ 2m 
2m 



(1) 



2 m 



k\ i\-2m) u 1 ; (§)„ V 2 2 

k=0 V2 Ik \2/2m x 



Rummer's summation formula (see [551 p. 10]) implies (later we set a = —2m, b = |' 



2/' 1 («■ b: 1 i r, - I ) = 2-' ,i\ ( |, - 6; 1 + a - 6; 1 
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Now we set a = —2m and use the terminating form of Gauss's sum (the Chu- 
Vandermonde formula) to obtain 

2 F 1 (-2m, b; 1 - 2m - b; -1) = 2 2m 2 F 1 (-m, -m + ^ - b; 1 - 2m - b; 1 

il-™) , (I) 



l-2m-6) TO (b + m) 1 

Thus, setting b = | we obtain, that the integral equals 

1 t n g 3 (t)dt- * ^ 2m ^ m i2m)] 



i i + 2m (l + m) m m!(2m + l)(2m + 3)...(4m + l)' 

We use the variance of a complex random variable to give quantitative insight into 
Figs. [2]j4j The calculations leading to the following formulae are based on the values of 
integrals of monomials over the real or complex unit spheres. Suppose A is a normal 
matrix, written in the form A = VDV^ where V is unitary and D is diagonal with 
entries Ai, A2, • • • , A at- Thus the columns of V are eigenvectors of A. Let X denote the 
random variable (-01 A\ip) and let m = j^trA = Ylf=i If ^ is a random vector in 
the unit sphere in M. N or then SX = m (where S denotes the expectation). Denote 
the complex variance v (X) =£(|X-m| 2 ), a measure of (2-dimensional) spread of X. 
For the complex shadow one finds 



N 

— ml 



v (x) = — - Via,- 

K J N (N+ 1 J 

j=i 

For the real shadow and when A can be diagonalized by a real orthogonal matrix, that 
is, V is orthogonal (VV T = I) then 

2 N 

v (X) = — -> |Ao — ml , 

V ; N (N + 2) ^ 1 J 1 ' 



which is larger than the variance of the complex shadow (see Fig. 3(b) ). In the general 

2 

real case V is not orthogonal, then set C (V)^ = J2k=i VkiVkj for 1 < i, j < N. Thus 
C (V) is a symmetric unistochastic matrix, and its eigenvalues all lie in [—1,1]. By 
straightforward computations we find 

« w = N{N+2) E ^ - m ) ^ - m ) (*« + c 

a positive quadratic form in (A, — rn)^ =l . The eigenvalues of I + C (V) lie in [0, 2], and 
all equal 2 when V is real orthogonal. 

As example consider the unitary Fourier matrix whose entries are primitive N th 
roots of unity : V jk = exp (j - 1) (k - 1)) , 1 < j, k < N. Then C (V) ifc = 1 
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if j + k = 2modiV otherwise C {V)j k = 0. The eigenvalues of / + C (V) are 2 with 
multiplicity [yj + 1 and with multiplicity |_^-^J- The matrix C (V) also contains 
information about approximating the eigenvectors of A by points on S N ~ 1- . indeed 

min [Ef=i ~ cVjk\ 2 : \c\ = 1, ^ e S"' 1 } = 2- ^2 (1 + C (V) 1 ^. In particular if 
some Afc is an extreme point in the convex hull of {Xj : 1 < j < N} then is in the 
real numerical range of A if and only if ^7=1 Vjk = 1- This is illustrated by Fig. 



4. Product numerical shadow 



Consider the shadow restricted to the set of pure product states. More formally, we 
assume that Hd — %N ® apply the definition Q and take 1Z = {\ipA) <8> \^Pb)}, 
where \ipA) £ Hn, while j^s) € Ha/ and these states are normalized. The set 
1Z of separable (product) pure states has the structure of the Cartesian product 
K = CP^ 1 x CP M ^ c <CF MN -\ 

The simplest case of iV = M = 2 corresponds to the two-qubit case. The set of 
separable (product) pure states has then a form of the Cartesian product of two spheres 
S = S 2 x S 2 . In other words this set forms a Segre embedding, CP 1 x CP 1 C CP 3 . 

One may also consider the shadow with respect to real separable states. In the 
two-qubit case this shadow corresponds to a projection of the product of real projective 
spaces HP 1 x HP 1 c HP 3 , which forms a torus S 1 x S 1 = T 2 . Such a structure can be 



recognized on some plots shown in Fig. 5 (d) - (f ) 



The matrices used to obtain projections are following 
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One can note that the numerical ranges presented in Fig. [5] have a particular 
structure. Matrix £> 4(l is permutation equivalent to a simple sum of two matrices so 
its numerical range forms a convex hull of two ellipses one with focal points {1,-1} and 
the other {i, — i}. As these ellipses do intersect their convex hull contains four interval 
segments. A similar situation occurs for the matrix B^ where two ellipses - one with 
focal points {!,— i} and the other {i, — 1} - do not intersect and the convex hull has 
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Figure 5. Product numerical shadow for illustrative operators of size N = 4 with 
respect to complex separable states, panels |(a)| |(b)| and (c) form projections of 
(DP 1 x (DP 1 C (DP 3 , while projections with respect to real separable states shown 



in panels (d) (e) and (f) present projections of RP 1 x RP 1 equivalent to the torus 
T 2 . Dashed circle represents the image of the outsphere while dotted line denotes the 
standard numerical range. Plots are done for matrices translated in such a way that 
their traces are equal to zero and suitably rescaled as described in [28] . 



only two flat lines. For the matrix the numerical range is a convex hull of an ellipse 
with focal points at eigenvalues {1, — i} and the line segment between two eigenvalues 
{i, —1}. A similar situation arises for the matrix B ic . 

Consider a particular case of an operator X with the tensor product structure, 
X — A <g) B. Then its product numerical range is equal to the Minkowski product of 
numerical ranges A® (A (g) B) — A(A)^A(B) - for more information see [21] . 

Is in this case the shadow of A <g> B restricted to product states can be expressed 
by the numerical shadows of both operators, 



»%b(E) = [ ( f J dfi A (t) + l {0eE}f i A ({0})- 

J<c\{o} V 1 J 



(17) 



Here ha and \xb denote the probability measures related to the numerical shadows of A 
and B respectively, while E C C is a measurable subset of the complex plane. 
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4-1. Mean and variance for the separable numerical shadow 

In the case of separable shadow of matrix X G M N 2 it is possible to obtain explicit 
expressions for the mean and the variance. We have 

E (z) = ^trX , (18) 

and 

E W = N*(N + iy (' tr X|2 + 11 tr ^ X|| HS + II tr * X Hhs + II^IIhs) • (19) 
Formulas involve the partial traces tr^X and trgAT, follows from more general fact given 



in Appendix Appendix B These expressions imply directly the following result. 



Var(z) = N2 ^ + i)2 (-^^1 tr X| 2 + || tr A X||^ s + || tr fl X\\ 2 W + \\Xf HS j. (20) 
4-2. Separable numerical shadow for a diagonal matrix 

Consider a diagonal matrix X defined on a composite Hilbert space, T~Ld — 7~Ln ® 

of dimension D = NM. Its diagonal elements forming the spectrum {xi\f =1 can be also 

represented by two indices, {y^} with /j, — 1, . . . , N and u — 1, . . . , M. 

Let |1,1) be an arbitrary fixed pure product state in so the set of random 
separable pure states can be obtained as = U\l, 1) = W ® V|l, 1), where W G U(AT) 
and V G U(M) are independent random unitary matrices distributed according to 
the Haar measure. Thus the expansion coefficients of the product state read 
(U n , U 12 , ...,U 1D ) = (WnVn, . . . , W lN V lM ). 

The separable numerical shadow of the diagonal operator X is defined as the density 
distribution of random numbers z := (ip\X\if)), where is a separable random state 
defined by random unitaries U and V. In this case one has 



D N M 

z := (1, i|(iy W)x(ww)ii, 1) = Z^i* 7 "! 2 = J2Y,y^\ w ^\ 2 \ v ^\ 2 = x ~ r ( 21 ) 

i=l A*=l v=\ 

where r is a real probability vector of size D = NM. It can be considered as a tensor 
product of two probability vectors p G A A r_ 1 and q G Am-i, since its components read 
Tfj, v = Pfj,qv with with \i — 1, . . . N and v — 1, . . . , M. 

Thus the separable numerical shadow of a diagonal operator can be considered as a 
projection of the Cartesian product of classical probability simplices, Ajv_i x Am-i- In 
the simplest case of .0 = 4 = 2x2 the Cartesian product of two intervals (1-simplices) 
forms a square, which lives inside the tetrahedron of the 4-dimensional probability 
vectors. 

As in the previous case we can distinguish two probability measures in the space of 
unitary matrices. They lead to 
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C) complex separable shadow, generated by the Haar measure on XJ(N) and U(M), for 
which both probability vectors p and q are distributed uniformly with respect to 
the Lebesgue measure on the simplices Aat_i and Am-i, respectively. 

D) Real separable shadow, generated by the Haar measure on the orthogonal groups 
O(N) and O(M), which lead to the statistical measure (Dirichlet measure with 
s = 1/2) in both simplices. 

Note that in these case the separable shadow of X is supported on its product 
numerical range |21| . which in general forms a proper subset of the convex hull of the 
spectrum. The product structure of the classical probability vector r in (21 ), generalized 
for a multiple tensor product structure, is consistent with the parametrization of the 
product numerical range described in Prop. 12 in [21] . 



5. Maximally entangled numerical shadow 

Consider an operator X acting on a Hilbert space with a tensor product structure, 
H = Ha ®Hb- For simplicity let us assume that the dimensions of both subspaces are 
equal to iV so the total dimension reads D = N 2 . Among all pure states of the N x N 
system one distinguishes the set £ of maximally entangled states. It contains the states 
equivalent with respect to a local unitary operation Ua ® Ub to the generalized Bell 
state, \ip + ) = ^= J2i \h Thus the set of maximally entangled states has the structure 
of U(iV)/U(l) = SXJ(N)/Z N , where Z N is the discrete permutation group [32], [HI Ch. 
15]. Choosing 8 for the set R in Q we define the shadow Px(z) of an operator X with 
respect to the maximally entangled states. The corresponding probability measure will 
be denoted as djj, x (z). 



5.1. Two qubit case: D = 2x2 



In the simplest case of 2 x 2 Hilbert space the set £ of maximally entangled states has the 
structure U(2)/U(l) = HP 3 — see [8]. Hence the numerical shadow of an operator A 
of order four with respect to the complex maximally entangled states can be considered 
as a projection of the real projective space on the plane — see the shadow for some 
illustrative operators presented in Fig. 6 



If one considers a further restriction and studies the shadow with respect to real 
maximally entangled states, the result can be interpreted as an image of the space 
0(2)/0(l) = HP 1 = S 1 Observe that the illustrative shadows obtained in this case and 



presented in Fig. 6 (d) - (f ) show indeed projections of a circle onto the complex plane. 



In the special case of a diagonal operator B of size four its shadow with respect 
to complex maximally entangled states can be identified with a standard shadow of 
a reduced operator B' of size 2. This fact is formulated in the following proposition, 
proved in Appendix Appendix A 



Proposition 3 Consider an diagonal matrix of order four, X = diag(di, d^, ds, d^) 
which acts on a composite Hilbert space H = Ha <8> Hb and the reduced matrix 
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Figure 6. Entangled numerical shadow for illustrative operators of size N = 4 with 



respect to complex maximally entangled states, panels (a) - (c) form projections of 
U(2)/U(l) ~ RP 3 , while projections with respect to real entangled states shown in 



panels (d) -(f) present projections of 0(2)/U(l) equivalent to a circle S 1 . Plots are 
done for matrices translated in such a way that their traces are equal to zero and 
suitably rescaled as described in [2"5j . 



Y = Y(X) = ^diag(c?i+(i4, c? 2 +c?3)- Then the numerical shadow Px( z ) of X with respect 
to complex maximally entangled states S is equal to the standard numerical shadow Py(z) 
of the reduced matrix Y(X) of order two. 



5.2. Two quNit case: D = N x N 

In this general case the set £ of maximally entangled states forms a manifold of (iV 2 — 1) 
real dimensions with the structure of U(iV)/U(l) [8]. In the case of the shadow 
with respect to real maximally entangled states the space under consideration reads 
0(iV)/0(l). 

To analyse the entangled shadow of a diagonal matrix we can perform initial steps 
used to prove proposition |3j It will be convenient to use a slightly different notation 
and consider a diagonal matrix X of size D = N 2 acting on a tensor product space 
W = Ha ®T-Lb with entries Xuj^a^ where i,j = 1,...,N. Consider a local unitary 
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matrix £g> U^ B ' where U^ A \ E XJ(N). The entangled unit state v has entries 



1 N 



(A)rr(B) 

ik u jk 



N 



U, 



(22) 



where U = (U^) . For purposes of numerical shadow the integration with respect 
to the Haar measure over both matrices and U^ B > can be replaced by a single 
integration over the N dimensional random unitary matrix U. Observe U \- > U is a real 
automorphism of U (N) invariant for the Haar measure. 

Taking a diagonal matrix X^j^^j^ = 5 i 2 ) 5 (ji, j 2 ) we find its 

expectation value for random pure state v 

N 



v^Av 



\Ua\ 



N 



tr (CB T ) 



Here we reshape the diagonal matrix X of order N 2 to get a matrix C of order N with 
entries Cy, while .B stands for a unistochastic matrix, 5^ := l^l 2 , for i,j = 1, . . . , N. 
The case N = 2 studied above relied on the simple nature of the set of unistochastic 
matrices of order two, equivalent to an interval. It is known that the structure of 
unistochastic matrices for iV > 3 is complicated and interesting J32J [33]. Thus we are 
not able to formulate a direct generalization of proposition [3] for the N x N problem. 
For the case of 3 x 3 matrices consider the function 

1 



v'Av 



tr (CB T ) 



where the (variable) unistochastic matrix reads 



B 



h 

h 

h - 63 



bi 

b 2 - 



- 

I-63- 
64 h + b 2 + 63 



b, 



Y2k=l Ciki 



Given the matrix C let C, 
C. = Y^H,j=i Cij- For simplification we set 

1 



C, 



Y?k=i C kj, for 1 < < 3, and 



C, 



3 Cj 



1 



a. . 



3 3 J 9 

One can observe, that the row and column sums of [7^] are zero and 



v^Av 



-tr 



Ah + 2b 2 + 26 3 + 64 - 3 26i + 46 2 + 63 + 26 4 - 3 
26i + 6 2 + 46 3 + 264-3 64 + 26 2 + 26 3 + 46 4 -3 





7n 


712 




721 


722 



+ 9 C - ' 



Without loss of generality we can assume C = 0. Powers of v'Av can be integrated 
using the formula from [311 Prop. 3.3]: 



f b?b. 

JV(3) 



^Kfbfb^dm 



nx\n 2 \n 3 \nA(2) n2+m (2) ni+n . 



2+^4 



(3) 



ni+n2+ri3+"4 



(2) 

ni+ri2 

(2) 



4-^3 



-ni, -n 2) -n 4 , 1 + n 3 
1, 2 + n 3 , -1 - m - n 2 - n 4 : 
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where dm denotes the Haar measure on U(3). 
For example 




There is an interesting special case when one of the variables 6j, say 64, does 
not appear in v^Av. The triple (6i,6 2 ,6 3 ) is a point in the pyramid with square 
base {(0,62,63) : < 62,63 < 1} and vertex (1,0,0). The induced measure (from U(3)) 
on the pyramid is jz^dbi db 2 db 3 . The numerical range of A is an affine image of 
the pyramid, hence is the convex hull of the images of the vertices of the pyramid, 
{(1,0,0), (0,1,0), (0,0,1), (0,1,1), (0,0,0)}, that is, a convex polygon. For three 
arbitrary complex numbers -21,-22, £3 let 



711 


712 




—6-2i — 6z 2 


3zi + 3,2 3 


721 


722 




3z 2 + 3z 3 


-3z 3 



then v^Av = (z 3 — 2z\ — 2z 2 ) 61 + (z 3 — z\) b 2 + (z 3 — z 2 ) 63 + Z\ + z 2 — z 3 . The vertices of 
the pyramid are mapped to {zi, z 2 , z 3 , —z\ — z 2 , z\ + z 2 — z 3 }. It is possible that these 
points form a (not regular) pentagon. For example let Z\ = 1, z 2 = i, z 3 = — | + | then 
the range is the convex hull of {l, \ + |i, i, — | + |, — 1 — i}. 

It may not be easy to find the shadow measure explicitly, but one expects a higher 
density in the neighbourhood of —z\ — z 2 , the image of 6 = (1,0,0). Fix e > and 
consider the set B £ = {(61, 6 2 , 63) : 1 — s < 61 < 1,0 < 6 2 , 63 < 1 — 61}. The normalized 
volume of B £ is e 3 and the U (3)-measure of B £ is e 2 , so the relative density is -. 



5. 3. Mean and variance for the entangled shadow 

It is possible to get explicit expressions for the mean and the variance of the entangled 
shadow d[i s x of matrix X acting on N x N Hilbert space % = Ha The following 

results 

E(z)= [ z dfi £ x (z) = ^trX , (23) 



c 



and 



= J ^ *A ( z ) = N2 (^2 _ 1} { Whs + I trX| 2 } 



- N , ( N2_ 1) {\\ ^X\\l s + || tr B X||^ s } , (24) 

which involve the partial traces tr^X and tr#X, are derived in Appendix Appendix B 
from more general fact. These expressions imply directly the following result. 
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Proposition 4 The expected squared distance from the mean with respect to the 
entangled shadow d\i x reads 



[ \z-E(z)\ 2 dp x (z) 
J<c 



1 



N 2 (N 2 - 1 
1 



I^IIhs + 



1 

N 2 



trX\ 



^{11^X11^ + || tr fl X|||s}. (25) 



N 3 (N 2 

Let us apply these formulae in the special case of a diagonal matrix Xf^j^u^j^ = 
8 (21,^2) 8 (ji, 32) Ch,h- In this case the necessary ingredients of (25) simplify 
considerably, e.g. 



N 



N 



i,j=l i,j=l 

N N 



•3 ' 



N 



tr((tr A X) (tr A X^))=J2 



i=l 



N 
i=l 



JV 



tr((tr s X) (tr B Xl))=J2 



i=l 



N 



i=i 



For the easy case N = 2 we find E [\z — E (z)\ 2 ) = ^ \Cu + C22 — C12 — C21 1 2 , in 
agreement with the previous calculations. 



6. Shadow and dynamics of quantum entanglement 

In previous sections we analysed the entire set of quantum states with its subsets and 
their projections onto a plane. In this section we specify a concrete quantum dynamics 
(in general non-unitary), choose an initial quantum state p(0), and following [28J we 
analyse its trajectory projected on the plane of a shadow of a selected non-Hermitian 
operator A. In particular we will be interested in dynamics of quantum entanglement, 
so the separable shadow of A will be used as a background for the trajectory obtained 
from the expectation values, z(t) = Tr Ap(t). 

Investigation of the dynamics of quantum entanglement was initiated in [10], in 
which the evolution of certain measures of entanglement in time was studied for a 
model non unitary dynamics of a two qubit system and several qualitatively different 
scenarios of behaviour of entanglement in time were identified. In particular, revivals 
of entanglement in time and an effect of sudden decay of quantum entanglement was 
reported. The latter effect was later called entanglement sudden death by Yu and Eberly 
[4T| |4"2"] and the dynamics of entanglement was studied by several authors in various 

setups [131112112 113 SI]. 

In general, all these dynamical effects can be explained in a simple geometric 
manner, if one takes into account the structure of the set of separable states analysed 
in [2]. As the convex set <2^ p of separable mixed states of a bipartite system occupies 
the central part of the set Qe> of all states of size D, entanglement revivals occur if the 
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Entanglement 
sudden death 



(a) Schematic figure often used to explain dynamics of quantum 
entanglement 





(b) Actual trajectory observed with the 
separable shadow of an exemplary matrix 
Xi in the background. All points outside 
this shadow are entangled (red crosses along 
the trajectory), while the points projected 
into the centre of the separable shadow are 
separable (blue circles along the trajectory). 



(c) The same trajectory observed by matrix 
X%. States of the trajectory projected into 
the separable shadow are typically separable. 



Figure 7. Dynamics of quantum entanglement. 



unitary dynamics moves the initial state several times across the separability boundary. 
On the other hand, entanglement sudden death effect takes place, if the decoherence 
is so strong that the initially pure state gets mixed in such a pace that it crosses the 
separability boundary only once. These possible scenarios are shown on a schematic 



sketch, used in a conference talks for several years — see Fig. 7(a' 



Making use of the technique of numerical shadow we are now in position to observe 
similar behaviour of entanglement for a concrete choice of quantum dynamics and initial 
states. Consider the following discrete-time dynamic quantum process of a system 
consisting of two qubits [10]. The system is assumed to be initially described by a 
maximally entangled pure state, p — |(|00) + |11))((00| + (11|). The discrete time 
evolution of the system is given by one-step unitary evolution expressed by Pauli 
matrices, U = e ia<Tx ® ay , followed by an action of the depolarising channel acting locally 
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on the second qubit, 



Hp) 



(3t 2 <g> cr p , 



(26) 



p£{x,y,z} 



so that pt+i = U$(p t )U ] . 

As the initial state is chosen to be pure, the trajectory begins at the boundary 
of the set of mixed states. As time t increases the trajectory plunges into the set of 
mixed states periodically crossing the set of entangled states. Artists impression of such 



processes is depicted in Fig. 7(a) where the outer oval corresponds to the set of the pure 



states, its interior corresponds to the set of all mixed states, ellipse corresponds to the 
set of separable states and spirals depict the trajectory. 

Let us fix the parameters of the discussed process, by setting the interaction strength 
a = 0.1 and decoherence rate /3 = 0.03. We arbitrarily chose two matrices 



X 1 



-1 i 

10 

0-10 

1 



and X 2 



1 i 
0-100 

0-10 

1 



(27) 



which allow us to project the trajectory p t onto the complex plane. We can calculate 
the images of the trajectory zj® = tr(p\Xk),k = 1,2 and superimpose them on the 



separable shadows of matrices X^. The resulting images are shown in Fig. [7| panels |(b) 
and (c) The red crosses indicate entangled states and blue circles indicate separable 
states. 

It can be easily seen that by choosing an appropriate observation matrix it is 
convenient to observe the dynamics of entanglement in the process. In the general 
case for any given trajectory, it is hardly possible to find such an observation matrix 
whose product numerical range contains only images of separable states, but images 
of entangled states always lie outside of product numerical range. Therefore product 
numerical range and separable shadow of a matrix are useful tools to visualize dynamics 
of entanglement, the effects of entanglement sudden death and entanglement revival. 



7. Multipartite systems 

It is natural to ask about the properties of the numerical shadow in the case when one 
aims to study composite quantum systems consisting of more than two subsystems. 

Let us consider the simplest case of a multipartite quantum system, i.e. a 
system composed of three qubits. In this case N = 2 3 = 8. As an example we 
take an unitary matrix Us of size eight written in the standard computational basis 
{|0,0,0),|0,0,1),...,|1,1,1)} 



2i7T 2i7T 2i7T 2i7T 2i7T 2i7T 

U 8 = diag(l, e3,es,e 3 ; e 3 , e 3 ( e 3 ; 



(28) 
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Figure 8. Restricted shadows given by matrix from Eq. 28 In each case the 
probability distribution is supported on a subset of a numerical range, namely restricted 
numerical range. In the general case this subset is not convex and the restricted 
numerical shadow can be supported on non convex set. 



The product numerical range of this operator is not simply connected [TjJl [21] . so 
it is instructive to study the shadow of U$ with respect to the space of 

a) all pure states (standard shadow); 

b) product states (product shadow), \ip) sep = Ua® Ub <8> Uc\0,0, 0); 

c) GHZ entangled states, |V)ghz = U A ®U B ® £7 C (|0, 0, 0) + |1, 1, l))/V2; 

d) so-called ^-entangled states, |^) w = U A ®U B ®U C (\1, 0, 0) + |0, 1, 0) + |0, 0, 1))/V3- 

Here Ua, Ub, Uc are independent random unitary matrices taken from U(2) with respect 
to the Haar measure. Fig. [8] presents the shadows generated by matrix U$ with respect 
to those classes of states. 

Observe that a generic operator acting on the three-qubit system leads to different 
shadows, if they are taken with respect to GHZ states and the W states. This is 
a consequence of the different topology of the orbits with respect to local unitary 
transformations produced by these two classes of entangled states. Study of numerical 
shadows restricted to certain classes of entangled states can contribute to a better 
understanding of the geometry of the manifold of locally equivalent states. 
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Investigations of the restricted shadows of operators acting on multipartite systems 
lead to a wide class of interesting problems which in general are difficult to solve. 
However, some results obtained in the previous sections for the bipartite setup can 
be generalized for multipartite systems. 

Consider an operator acting on the composite Hilbert space describing m-partite 
system of dimensions N\, N 2 , . . . , N m , respectively. Its shadow with respect to separable 
states |^) = <S> \4>2) ® ■ ■ ■ \ipm) leads to classical product measures on a simplex of 
a composed dimension D = N1N2 ■ ■ ■ N m — 1, induced by the Dirichlet measures on m 
simplices of the dimension (iVj — 1) with i = 1, . . . ,m. 

8. Concluding remarks 

In this work we propose to combine the notion of restricted numerical range with the 
numerical shadow of an operator which is a probability measure on the complex plane. 
On one hand the numerical shadow can be investigated for a given operator X. On the 
other hand, one may analyse the shadows of all operators of a fixed dimension D, which 
can be considered as projections of the set of all pure states of size D onto a plane. 

In a similar way it could be interesting to study the shadow of a given operator 
X with respect to various sets of pure states. For instance, in this work, we analysed 
the standard shadow with respect to complex states and the shadow restricted to real 
states. For operators acting on composite Hilbert spaces one can additionally study the 
shadow with respect to separable or maximally entangled states, complex or real. Note 
that in general the restricted numerical range is not convex f2M |2"5] , which implies that 
the restricted numerical shadow can be supported on non convex sets. 

Following the complementary strategy one may take the set of all operators of 
a given dimension and analyse their shadows restricted to a certain class of states. 
These probability distributions on the complex plane convey some information about 
the structure of these specific subsets of the set of all quantum states. Consider for 
instance the simplest Hilbert space with a tensor product structure Hi ® Hi- Then the 
standard shadow carries information about the set of all pure states which forms the 
complex projective space CP 3 , while the shadow restricted to real states corresponds 
to the real projective space HP 3 . 

In a similar way, the shadows of matrices of order four with respect to separable 
states illustrate the projections of the product of two spheres CP 1 x CP 1 = S 2 x S 2 , 
while the real separable shadow corresponds to projections of a torus RP 1 x HP 1 = 
S 1 x S 1 = T 2 . The shadow with respect to maximally entangled states (also called 
briefly the entangled shadow) represents the set £ = U(2)/U(l) also equivalent to real 
projective space RP 3 . In the case of real entangled shadow we observe projections of 
the set £ R = 0(2)/0(l) equivalent to the circle RP 1 ~ S 1 . 

The notion of separable shadow is useful to analyse the dynamics of quantum 
entanglement. For a given initial quantum state p(0) and a certain dynamics we select a 
non-Hermitian matrix A and study trajectories on the complex plane formed by the time 
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evolution of its expectation value, z(t) = Tr Ap(t). Investigations of such a trajectory 
superimposed on the separable shadow of A contribute to our understanding of the 
dynamics of quantum entanglement and allow us to visualize the effects of entanglement 
sudden death and entanglement revivals. 

The notion of restricted numerical range is easily formulated for operators acting 
on Hilbert space with multiple tensor product structure [25J, which correspond to 
multipartite quantum systems. Therefore it is natural to define restricted numerical 
shadow for various classes of quantum states of multipartite systems. For instance, in 
the simplest case of a three-qubit system, described in the Hilbert space H,g — H.2 
one may distinguish two classes of maximally entangled states called GHZ and W, 
which cannot be locally converted in any direction [35J. Studying numerical shadows 
of matrices of order 8, restricted to GHZ states or W states, we are thus in position 
to investigate the differences between the structure of these two important classes of 
three-qubit quantum entangled states. 

In conclusion, we have introduced the notion of the restricted numerical shadow 
of an operator and established its basic properties. On one hand we believe that this 
topic is interesting from the mathematical point of view, as it relates operator theory 
and probability. Moreover, we are tempted to expect that further investigations of the 
restricted numerical shadow will contribute to a better understanding of the geometry 
of quantum entanglement, so that they become directly applicable to the theory of 
quantum information. 
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Appendix A. Proof of Proposition [3] 

To analyse properties of the shadow with respect to complex maximally entangled states 
for an operator acting on the 2x2 Hilbert space let us analyse the structure of a local 
unitary transformation <8> acting on H = Ha ®Hb- Consider two generic 
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e i(V>i+0i) cos Q 1 _ e i(V>i-02) sm Q x 

e Wi+fe) sin 1 e Wi-4>i) cos e 1 



gbh+te) cos 6 2 sin 9 2 

e i(^+«4) s j n ^ e i(^2-^>3) cos # 2 



where — tt < <pj < n, < ^ < n, < #j < | and the Haar measure (for C/^) is 

dm(U) = -—rdipid(j)id(j)2 sin #1 cos9\d9\. 

Taking tensor products of corresponding columns of and £/ < -' B - ) and adding them we 
obtain a parametrization of an entangled state, 

tt( a )tt( b ) , tM)tt^ 



V2 



Consider an diagonal matrix of order four, X = diag(di, d 2 , d 3 , d 4 ) 
so that its expectation value for a maximally entangled state reads 

v ] Xv = ^ (di + rf 4 ) 9i + ^ (d 2 + d 3 ) g 2 , 

91 = cos 2 #1 cos 2 6 2 + sin 2 sin 2 6 2 H — sin 2^ x sin 2# 2 cos (0i + 2 + 03 + 04 ) , 

2 

g 2 = cos 2 61 sin 2 # 2 + sin 2 9\ cos 2 # 2 — ^ sin 2#i sin 2# 2 cos (0i + 02 + 03 + 04) • 

Observe that q\ + q 2 — 1. Compare this situation to the complex shadow of the reduced 
operator Y = Y(X) 



Y 



\ (di + d 4 ) 

§ (d 2 + d 3 ) 



x^Yx = \ (di + d 4 ) |xi| 2 + ^ (d 2 + d 3 ) |ic 2 1 2 , 



where x belongs to the unit sphere S 1 in C 2 . We know that 

n\k\ 



L 



I 1 2n 1 \2k 1 / \ 
\X\\ \X2\ djJ,[X) 



(k + n + 1)! 

To identify the two shadows it only remains to show 
q n iq k 2dm{U^) dm{U^) = — nUA 

'U(2)xU(2) 



, k, n — 0, 1, 2, 



-, fc,n = 0, 1,2, .... 



(n + fc + 1)!' 

In fact it suffices to show f g™ = ^j-j-. The integration over the angles 0j can be combined 
into one 0, due to the rotational invariance. Write 

qi = (cos 0i cos 9 2 + e l<?i sin sin 9 2 ) (cos #1 cos 9 2 + e~ l * sin sin 9 2 ) , 



i=o 1=0 



n \ in 



(cos 0! cos 9 2 ) 2n ~ ] ~ l (sin 0! sin 9 2 ) j+l e l{j - l)(t> . 
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Now integrate with ^-d(j) over — tt < <f) < it, then over the dd\dd2 part: 

j q l = 4 (j) J 2 J 2 (cos #i cos # 2 ) 2n - 2j+1 (sin B x sin # 2 ) 2j+1 dM# 2 

fMV ("-j)!j! ) 2 _ 1 Vx- 1 



v /7 V (n + i)l 7 (n • lrf^ « + 

To complete the proof 

k 



j=0 

h 



V (k\ ( _ iy 1 = _J_ (zgj (" + ^ 
n + j + l n+lj^ j\(n + 2) j 

(n + 2-n-l) A n!A;! 



(n + 1) (n + 2) fc (n + fc + l)!' 
where is Pochhammer symbol. This used the Chu-Vandermonde sum (terminating 
2 F l(1 )) a „d 1 - + - {n+l) ' 



n + l+j (n + l) j (n + l+ j) (n+l)(n + 2) j 
Appendix B. Expectation values for shadows 

Consider a general matrix X acting on a iV x N composite Hilbert space H = Ha^Hb 
with entries written in a four index notation , where ji, i 2 , ji — 1, ■ ■ ■ ,N. Upper 

pair of indices determines the row of the matrix, while the lower pair determines its 
column. We use the standard operations on matrices, which in this notation read 

TV 

X} ii2 :=X jlh , trX:=J2 X v ( R1 ) 

and introduce reduced matrices of size N obtained by a partial trace over a single 
subsystem, 

N N 

(tr A X) nj2 := ^I ii21 (tr B X) nn := ^l 8lJ (B.2) 
i=i ,=i 

By the formula in Collins and Sniady [34J we have 

UijUi'fdm (u) = -rfSn'Sjji, (B.3) 

U(JV) ^ v 

Ui^Ui^Ui' j'Ui' f dm (u) = 1 {Shi'^i' Sj^S^f + Si^ Si^Sj^S^f } (B.4) 

U(7V) iV — 1 

~~ N (N 2 - 1) {^'1^2 'Wa'W + ^ 2 ^AAiJ > 



Restricted numerical shadow and geometry of quantum entanglement 



29 



where dm denotes the Haar measure on U (N). 
Let 

z = ( X \U® VXU ] ® V ] \x) 



(B.5) 



where U, V are stochastically independent random unitary matrices of size N, 
distributed with Haar measure and \x) is an arbitrary vector. We are interested in 
obtaining mean and variance of z, thus we will calculate E(z) and E(zz). We have 



E(z) = E((x\U g> VXtf <8> V ] \x)) = (x\E(U ® VXE/t <8> 



(B.6) 



and 



E(zz) = E({x\ <8> <x|(Z7 <8> VXf/t <8> F f ) ® (17 <g> VX[/ T <8> V T )|x) <8> (B.7) 
= (ac| <g> (z| E((C/ ® VXE/t (8) V*) ® (C7 (8) VXU T ® \/ T )) |ar) ® |x)(B.8) 
To obtain the mean we calculate 



.W=l 

AT 

/ E(u il fc 1 f j 2 fc 2 Xfc 1 fc 2 ; Uj 1 i 1 



^2, 



(B.9) 
(B.10) 



, ki,k2,h,h=l 
N 



E 



jjjffiiihSkih8i 2 j 2 h2i 2 Xk 1 k2 } < B. I I ; 



AT2 



trX 1 



N 2 xN 2 > 



(B.12) 



where we have used a convention, that we split indices £ = iV(£i — 1) + £2, and £1,^2 
has values from {1, 2, . . . N}. From the above, we have that for any normalized \x) 



E(z) =E((x\U ®VXU ] ® V ] \x)) = ^trX. 

N 2 

Let us now we calculate the second moment. Let 

M = E((U (8) VXU ] <8> V ] ) ®{U® VXU T <8> V T )) 



We have 



(ii,i 2 ,A;i, /c 2 |M|ja, j 2 ,h,k) 



iV 2 



/ . E(tt ilQ , 1 t'j 2Q , 2 Xaia2Uj l/ 3 1 Uj 2/ 3 2 Mfc 17l Wfc 272 X7l72Mz iei t'/ 2 



(B.13) 

(B.14) 

(B.15) 
(B.16) 



a,/3,7,e=l 
A' 2 



V" E(u iiai M /iei u il/3l u fcl7 jE(t; i2a2 ^ 2e2 w i2( g 2 w A . 272 )Xaia2X7i72, (B.17) 

*— ' P\P2 e l € 2 

a,/3,7,e=l 



Using formula ( |B.4[ ) for the integral we have that M = c x M^ +c 2 M^ +c 3 M ( V T c 4 I (4) , 
where 
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N 

M«= \h,i2,k 1 ,k 2 )(ii,h,k 1 ,k 2 \, (B.18) 

ii,i 2 ,fci,fc 2 =l 

Af2 



Similarly 



1 



Cl 



c 2 



Next we have 



(iV 2 - l) 2 



^ (<WAi7i ~~ T^ai7i<Wi ) 



^a 2 #A 2 72 - ^a 2 72<W 2 ^ ^g^l}^ (B.19) 

I ( N - 1 N — 

, , TO —7: / v\"aict2^7l72 — 7 JTaia 2 JT7ia 2 

(TV 2 - l) 2 \ ^ Q l"2 7172 TV ^ c*lft 7l/?2 

\ ai, "2,71)72=1 ai,a2,ft,7i=l 



— — > -Xaic*2 ATai72 H — 777: 7 .Xaact2 AT«i(*2 

TV ^ fta 2 ft 72 A^ 2 ^ ftft /Sift 

ai,a2,Pi,72=l ai,o:2,Pi,ft=l > 

(iV 2 _ 1)2 1 tr(X)tr(Xt) - I tr(tr A (X) tr A (X)t) - I tr(tr B (X) tr B (X)t) 
+ ^tr(XXt)| (B.20) 
( I trX| 2 - 1|| Mollis - 1|| trellis + ^||X||| S ) . 



iV 



M(2) = \ji^2,ji,k 2 )(ii,i2,ii,k 2 \, (B.21) 

ii,»2 ,ji,fc 2 =l 

AT2 



fjV 2 — l) 2 ^ (^i7i^ftei jy^Qi/3i^l7i ) 

^a 2 ft^£ 2 72 - ^^a 2 , 72 ^ft,e 2 ^ ^jgX|i72 (B.22) 

' 'll Mollis - ^II^IIhs " 4l «| 2 + 4ll M^lllsl . 



(AT2 _ 1)2 ^11 /N«!> TV " Hb AT 1 A^ 2 

N 

M (3) = ^ \ii,j2,h,j 2 )(i 1 ,i 2 ,k 1 ,i 2 \, (B.23) 

ii,i 2 J2,fci=l 

^3 = T^^p (ll Mollis - iyXUls - 1| trX| 2 + -^|| tr B (X)||| s 
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Finally 

N 

M (4) = \ji,h,ji,j2)(ii,i2,ii,i2l (B.24) 

*l,*2,ii,j2 = l 

C4 = jj^y 2 (|| tr B (X)f m - i|| tr A (X)\\l s - i|| tr B (X)|| 2 s + ^| trX| 2 ) . 

The integrals of the above type, for fixed dimension, can be calculated with the use 
of computer algebra program IntU [37] . 



Appendix C. Variance values for shadows with fixed Schmidt numbers 

Let us consider any pure state \x) from N x N composite Hilbert space H = Wa ® Hb 
with fixed Schmidt numbers A 1; A 2 , • • • , \n, thus we have 

N 

\x) = 

i=l 



>=x;v^i <a >® i«' b >> (c.i) 



for some orthogonal bases \i A ), \i B )- For simplicity we can take computational bases. 
We have the following identities 

= 1 (C2) 

N 

= £a 2 . ( c - 3 ) 

1=1 



(x\ <g) 


(x|Mi|x) £ 


?> x) 


= (x| <8 (x\M A \x) <8 | 


(x\ <g) 


(x|M 2 |x) £ 


3 x) 


= (x| <8> (x|M 3 |x) <8> | 


have, 








(x\ <g) 


(x|M|x) (g 


x) 


= Ci + c 4 + (c 2 + c 3 ) ] 



AT 

£A 2 , (C4) 

i=l 

Note that, the above depends only on a purity of reduced state. In the case of maximally 
entangled state (Aj = -k for i = 1, 2, . . . N) the second moment is given by 

(x\ ® (x|M|x) ® |x> = ^ _ (| trX| 2 + HXHy (C.5) 

1 



^3(^2 _ I) (H^P0IIhs + ||tr B (X)||| s ) (C.6) 



The above implies the result (25) for the variance with respect to the entangled shadow. 
In the case of separable states (Aj = 5i ti for i = 1, 2, . . . N) we obtain 

(x\ (8 (x\M\x) (8 |x) = (C.7) 

! " trX| 2 + || tr A (X)|| 2 HS + || tr B (X)||| s + ||X|| 2 S ) . (C.8) 



n 2 (n + iy 



Which implies the result (20) for the variance with respect to the separable shadow. 
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